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Abstract 

Closed form solutions are presented for sound 
propagation from a lino source In or near a Bhoar 
layor, The analysis is exact for all frequencies 
and is developed assuming a linear voloctty profile 
in the shear layer. Tilts assumption allows the 
solution to be expressed In terms of parabolic cy- 
linder functions. Tho solution La presented for a 
lino monopole source first embedded In the uniform 
flow and then In the shear layer. Solutions are 
also discussed for cartdin types of dipole and 
quadrupolo sources. Asymptotic expansions of the 
exact solutions for small and large values of 
Strouhal number give expressions which correspond 
to solutions previously obtained for tliCBO limiting 
esses , 


Introduction 


Bslsafl gives expressions for pressure fields of 
various high frequency convoctcd sources embedded 
in a unidirectional sheared flow*. The expressions 
include rhe simultaneous affects of fluid and source 
convection and refraction. Those results are for 
axlsymmotrlc flows with arbitrary velocity profiles . 

Tho next step In determining the effect of a 
shear layer upon round propagation Is to obtain 
analytic solutions for the far field pressure for 
the entire acoustic frequency range, ThlB has been 
accomplished In the present investlgatl)n^ at NASA 
Lewis Research Center by using an exact nnniytlcol 
solution obtained by Goldstein and Rice 1 * 5 not only 
for a monopolo source in a uniform flow near a 
llnenr shear layer but for dipole and quadrupote 
sources as well. This analysis Is further extended 
to the case In which tho source la embedded In the 
sheur layer for a monopolo, and a quadrupolo type 
of source for the entire frequency spectrum. 


Sound waved passing through a shear layer arc 
altered both in amplitude and direction. Hence an 
observer situated in the far field on the opposite 
« .<lc of a shear layer from a sound source will, in 
goncral, hear sound which differs from that origi- 
nally produced by the source. Any sound produced 
by an aircraft cnglno must p^nn through a shear 
layer or combination of shear layers. Consequently 
a thorough understanding of this phenomenon Is es- 
sential for determining the far field sound pres- 
sure amplitude directivity generated by an aircraft 
engine In flight. Such knowledge is also necesnary 
for developing adequate correction procedures for 
data from opcn-Jct anechotc wind tunnels used for 
simulating flight effecta on noise sources. 

There have been several studios of the effect 
of relative motion between two media on the propa- 
gation of sound waves. Rlbnor* ond Miles'* Inde- 
pendently Investigated plane waves Impinging on a 
velocity Interface between two moving streams. 
Gottlieb^ studied a sound source near a velocity 
discontinuity. Ho obtained far field solutions by 
evaluating the exact sound field integrals by the 
method of stationary phase, Slutsky, Tamagno and 
Horctti^- 6 extended Gottlieb’s analysis to include 
quadrupolo sources embedded In the menn flow. 

Their analyses wore developed for the sources situ- 
ated in nn axlsymmotrlc slug flow. Graham and Gra- 
ham' treated the propagation ot plane waves through 
a finite alipar Layer. The Bame authors 8 later in- 
vestigated the effect of a finite shear layer on 
sound waves from a point source. They obtained 
velocity potential solutions by means ol a aeries 
expansion about a singular point and presented re- 
sults for very low frequencies with acoustic wave- 
lengths much greater than the shear layer thickness, 

Goldstein^**® obtained solutions for low fre- 
quency sound from multipolo aourcca In axtsymmcLrtc 
shear flows with arbitrary velocity profiles. The 
results of his analysis show how tile mean flow af- 
fects the radiation pattern from the sources. 
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SymliolB 

A{k) arbitrary coefficient In equation Ci2) 

U(k) arbitrary coefficient In equation (14) 

b transformation parameter defined by equa- 

tion (IS) 

C,(k) arbitrary coefficients in equation (20) 

c q speed of sound, 340.46 m/scc 

f, arbitrary unit source term in uniform flow 

•* region 

f T j Fourier transform of f^ 

f externally applied force 

gj arbitrary unit source term in shear layer 

g^j Fourier transform of 

h coordinate of source location In -y di- 

rection, meters 



k, wave number In x direction In region 1, 

ko cos " 

k wave number In x direction 

H Mach number In shear layer 

Mach number ol uniform flow 
o 

F(y,k) Fourier tratialorm of p(y,x) 
p perturbation pressure, N/m^ 

q external volume flow source, equation (42) 

r radial coordinate of observation point In 

cylindrical coordinate system, meters 

b/ C 

S — ij — — Strouhal number (reduced frequency) 


ERRATA 

Figures 3 to 6 and 7(c). The dashed lines in the keys should be 

le versed, that is, the short-dash line denotes the high-frequency 
appi oximation, and the long-dash line denotes the low-frequency 
approximation. 


t time , aoc 

U velocity in x direction in the shear 

layer, moters/sec 


U 


Velocity of uniform flow, mctcrs/aoc 


V(y,k) Fourier transform of v(y,x) 

v perturbation velocity in y diroction 

x coordinate in diroction of uniform flow 

y coordinate 1 to direction of uniform 

flow 

2 defined by equations (13) and (IS) 

fl(y) dlrac delta function 

6 ahoar layor thickness, motors 

0 polar coordinate (polar angle) or direc- 

tion botwoon lino connecting source and 
obsorvotion points and direction of uni- 
form flow 

critical angle below which Bound waves 
will not propagate through shear layer 

/i wavelength 

£ transformation variable dofinod by equa- 

tion (1C) 

a 

P density 1.23 kg/m 

to angular frequency 

Subscripts: 
d dipole 


in 

1 


monopolo 

quadrupole 


2. In the uniform flow region U • U Q (region 2) 


2 ±i f 

P2 ’ c 2 jjj. 2 *2 ‘ f J 
0 


J - 1,2 (2) 


where 


a. u ^ 

Dt J TJt + U o Sx 


and 


f 6 <x)6 (y + h)e 

v{„ 


lu)t 


i - 1 

j = 2 


accounts for tho presence (J = 1) or absonce 
(j *> 2) of the source (of unit strength) is the 
uniform flow. 



D ,.2 X D ,A tTn 

n " ' ? ' c 2 Dt 2 ) ‘ 


J “ 1.2 
(3) 


where the prime denotes differentiation with re- 
speut to y, 


t|i a 


dll 

dy 


± -A +U ji 

Dt -it T u dx 


Analysis 


and 


A typical turbofsn engine is represented in 
Fig. 1. This indicates haw sound generated by such 
an engine might encounter vorious shear layers. 

A simplified formulation Of this problem La 
based on the flow geometry shown in Fig. 2 whore a 
medium moving with a unifonn velocity U (region 
2) and a medium at rest (region 1) era separated 
by a single shear layer having a linear velocity 
profile. Hie no flow case is used on ono aide of 
the shear layer for simplicity, and the solutions 
can he easily extended to the more general case of 
finite uniform flow on both sides. Cases will be 
considered in which a line source is located in tho 
uniform flow and In tho shear layer. 


(jqr(b006(y + hjo 1 ' 116 ) 


J “ 1 
J - 2 


accounts for tho absence (j - 1) or presence 
(J - 2) of Clio source In tho ahonr layer. 


The radiatlonship between tho transverse par- 
ticle velocity and tho acoustic pressuro Is given 
by 


Dt > 


Uy 


Tile governing equotions for sound propagation 
in each of these regions may bo found in Ref. 14 
(p. 10) and are given in terms of acoustic pres- 
sure as follows: 

i. In the no flow region, U *. 0 (region l) 


v 2 p, 


2 

" P» 


(1) 


The linear velocity profile in the shear 
layer is given by 


U(y) -U 0 y/6 

For the boundary conditions we hove the 
ocoustic radiation condition which requires out- 
ward propagating waves at ye ♦*., At the edges 
of tho sheor layer we require continuity of pres- 
sure 


p b pj at y » 0 and p * p, ot y -f. 


(4) 


T 


am) continuity of particle displacement which can 
be expressed In terms of the transverse perturba- 
tion velocity (slncu there Is no slip) 


■ Vj nt y *. 0 and v “ at y • -6 


(5) 


Assuming a simple harmonic time dependence and 
applying the Fourier transforms 


Jiwfc f " ikx 

p(y»h) * ~2v~ J cix 

•ou 

-lest r a 

,w - Vy v < x -y> 11 d * 

tfi 


v(y 


2,j •« ° lkx d * - i 


glvis the following system of second ordor ordinary 
differential equations and boundary conditions. 


In region 1 


p£ 4 


In region 2 


4-k 2 

L°o 


P : » 0 


4 


wliero 


ft ’ "» k ) ! 


P B f 

2 r TJ 


( 6 ) 

J “ 1,2 
( 7 ) 


12? 6 <* + 1,) 


Tj 


.! 1 
J » 2 


and 


H » U / c 
o oo 


In the sheer layer 
2H»k 


1 >" 




P‘ 4 


HH 


P “ R. 


TJ 


^ | 1 


2? 1 (y + W 


J = 1,2 (8) 

J - l 

J - 2 


where 



P' - 

I f 0 * 

-‘v. - 

kM^V 

(9) 

and 


H , JL 





c o 



with 

the boundary 

condition* 




P - P t 

J V « Vj at 

y » 0 

(10) 


P » P 2 

{ V = V 2 st 

y - -A 

01) 

The aolution of thla ayatem will dupond 
whether the source is located in tho uniform 

on 

flow 


or In tire shear layer. However, the procedure la 
the same in cither case. 

The solution of equation (6) in the no flow 
region which satisfies the radiation condition has 
the form 


A(k)u 


-v 


where 


h"V 


k 2 - (u)/c 0 ) 2 


( 12 ) 


(13) 


Wien the source La In Vho uniform flow region 
J » li f. » ferf>(y 4 h), g™ a 0 the solution to 
l l V r l 

equation (7) which satisfies the radiation condi- 
tion has tho form 


H(k)o 


V 




-y 2 [y+ii [ 


where 


72 




k 2 - (uVc o - kH o ) 2 


(M) 


(15) 


Tlie assumption of o linear velocity profile 
in the shear permits use of the solution obtained 
by Colds tcin and Rlcel 3 which exproaaes the acous- 
tic pressure in the shear layer in terms of the 
parabolic cylinder functions of Weber. 15, lb Hence, 
following Goldstein and Rlcel 3 , the new Independent 
variable 1 and dependent variable 11 are defined 
by 


*®ft'“) 


(16) 


and 


4 (•-**) «> 


Wliero 


2iH' 


(18) 
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H' 


These tranaformatlon* ara introduced into 
equation (9) which is then integrated oncu as in 
Ref. 13 to obtain Weber's equation. 


Putting the appropriate expressions for y, and 
for P^(li) and p' 1 ^) in the numerator of 
equation (23) yields , 


U" - <1/4 t 2 + b)U « 0 (19) 

where the primes now denote dlf forontation with 
reapoct to £. This equation has as its solutions 
an arbitrary linear combination of parabolic cy- 
linder functions U(b, * £), (Abramowlts 16 ), lhus 
the solution to equation (8) la written as an arbi- 
trary linear combination 

P ° C h (k)P + {i) + CJk)P~(0 
of the functions 

P (f) . .(»:( »■. .it)) 

* i . t. 


( 20 ) 

( 21 ) 


TTM 0 k (b + l/2)P<b - i/2) 


g'r + - cV 


(26) 


Tills expression for A(k) can be put into 
equation (12) for Pi which is tho Fourier Trans- 
form solution for the far field pressure pattern In 
region l. The actual pressure is then obtained 
from the Fourier inversion integral 


H* f -(ikx+y.y) 

Pl (x,y) o o / A(k)o dk (27) 


with A(k) given by aquation (26), Tlie above in- 
tegral is of the farm 


or 

P (?) -U(b - 1, 0 + (b - 1/2)11 (b + 1,-5) (22) 




dk 


Having obtained these solutions (Eqs. (12), (14) 
snd (20)) the boundary conditions (Eqs. (10) and 
(1.1)) are applied at the edges of the shear layer 
to evaluate the following expression for the coef- 
ficient A(k) in the far field pressure expression 

(Eq. (12)), 


A(k) a 


V 2 (-6+h) 


[p + (>, 1 )p , '((. l ) - p |H '(e 1 ?p'( f . 1 )] 

G'R + - G + lf 


wliere 

f(k) o lkx + Y^ (28) 

with Yi defined by equation (13). Tills typo of 
Integral can bo evaluated by the method of station- 
ary phase or sa ddle poin t method^! 18 (for large 

positive r ° + y 2 ) 


where 


o' = P' « 2 ) - y/(? 2 > 
R* * P' '(fj) + Y t P (^) 


(23) 


(24) 

(25) 


Tlie details of the evaluation of equation (27) 
are presented in Ref. 12. The result, given here 
in cylindrical coord liiaceu , is the following ex- 
pression for the far field pressure directivity for 
a lino monopolo source located in the uniform flow 


^Im 




2 6 


rM can 

ti 


»lll o 

(b + l/2)J'(b - 1/2) 


*- '. f v.'-vr-f- 

G R - R G 


and 


P , n y e .<»»!■ it). 

* l l J d‘ dy 


f = 1,2 (25) 


ti t , 


d£. 

dy 


■= -r21MI‘ 


where 


and 


" VK(0 


at y - 0 


1® « i 




2 2 

cos t ) -con v(-t.+h)+r+c t) (29) 




where 


■ (* € 


|(1 - M (J con >) + ^ (1 - H o cos n)^ - con^i'^ 

■ ll(b - 1, i, 2 ) 


+ /(l - H. cos • ) - )! (1 - M o cos u)- - con^f'^ 


(a - h 0 


4 


(36) 


x (b - l/2)U(b + 1,*^) (30) 

and 



k h 
K 6 


»•) 





- (1 


oln G)U(b 


l.Ht) 


+ (1 + aln fJ) (b - 1/2)1) (b + 1,H l ) (31) 


The frequency la defined in terms of Strouhal 
number ao follows : 

Recall that 



It is noted that when the source Is placed at 
the uniform flow odge of the shear layer, tn which 
case h “ 6 and £2 “ 4h» aquation (33) becomes 
equal to equation (29) • Both equations (29) and 
(33) rcduca to tho same special result regardless 
of source location when the line between the source 
and obsorvor Is oriented at 90° to tho Clow direc- 
tion (l.o,, 0 *• 90° In Eqs. (29) and (33)), This 

result is 


lm 



- ~ 


i 


8rri f- r 

<)«3Q° 

L C o J 


1/2 


-i^-(li- 26 +r)+lcet 


(35) 


This also agrees with Gottlieb's result^ at 

11 a 90° whan tho shear thickness A is act equal 

to zero. 


Approximate Solutions for l.tmltlnn Cages 


where k ° k^ 


k cos 11 was uacd. 
0 


of Lou and llli’li Frequency 


Note that h(on/c o )/ M 0 Is a nondlmonslonal group of 

parameters that occurs sevoral places in the pres- 
sure expression, llcncc wa define 


Strouhal number » S 



(32) 


Then 


and 


b ■» 1/2 1 S cos v 


f 


J: 


21s 


cos 1 • 


(1 - M coa 0) 


This Strouhal number is essentially a comparison of 
the shear layer thickness to the acoustic wave- 
length, and is referred to as a nondlmenslonal fre- 
quency. 


When considering the far field pressure for 
limiting values of frequoncy it is actutiiy limit- 
ing values of nondime no tons l frequency or Strouhal 
number (defined by Eq. (32)) which are of interest. 
Careful examination of the far field pressure equa- 
tions (Eqs. (29) and (33)) along with the defini- 
tions of b and £ reveals the strong dependence 
of the pressure upon StrouhoL number. This depend- 
ence tomes about through the parabolic cylinder 
functions which can be evaluated for small or large 
arguments by means of appropriate asymptotic ex- 
pansions found in Refs. 10 and 19. Tho parabolic 
cylinder functions appear only through A(k) as 
evaluated for k » k^ *» w/c a cos n. By evaluating 
these expressions for the limiting cases of very 
small and very large values of the Strouhal number 
the approximate expressions for the far field pres- 
sure will be obtained. These approximate exprea- 
stona for the far field pressure can also bo ob- 
tained by solving the differential equations asymp- 
totically for arbitrary velocity profiles ualng 
methods discussed in see 6,7,1 of Ref. 14. 


a 


Now placing tho source tn the slieor layer, 
2, f T = 0, e t ^ ~ t (y+h) In Eqs, (7) and 


1 2 " ’ “ L 2 2n 
(8)) and following exactly the same procedure 
yields the far field expression Cor acoustic pres- 
sure 


Pirn 


2s — r 
c 

o 


1/2 


(-1 sin >0 


t. h „ \2 

\ l ' 6 M o CD3 ) 


By replacing the parabolic cylinder functions 
in equations (29) and (33) with their asymptotic 
expansions* 6 far small arguments the low frequency 
approximations fur far field pressure are obtained. 
?or the source in the Uniterm flow 


lm 


2rr — r 
c 

o 


1/2 



- p eye + 

1 -i—r+iit 


(~i sin-*)e 0 


(33) 

r 2 

\ C R 

- g’ + r" J 

f- 

<U - M qpu ) 8ln« -r 


where 


COS L*) 2 'COB^l(-6+h)+t-C o t| 
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Far the aourco In the shear layer region 


Pirn * 


-il/2 


_ u> 

2ir — r 
0 


(l flin to (1 * M a coa 0 ) 

( i ■ * k o cds °y 


* < 


- 1— n list 
. C o 


- H q cob I>) 2 Dim' + f/(l - H q cob D) 2 - cos 2 0 


(37) 


bo □ f interest here, Tims far field pressure due 
to a lateral llnu dipole aourco in the uniform 
flow is the first derivative with respect to trans- 
vorae source location (h) of the raonopolu expres- 
sion (Eq, (29)). 

p ld ‘ "5T " ' r p lm ^ eo, *« * <1 - H fl cos 0) 2 

o 


(40) 

The lateral llnu quadrupole express ion la ob- 
tained by differentiating the monopolu expression 
(Eq. (29)) twice with respect to transverse source 
location (li) . 


It Should be pointed out that when the source 
ia located at the edge of the shear layer joining 
the uniform flow (6 = It) equations (36) and (37) 
reduce to the same reault. 



(1 


M coa u) 
o 


Replacing the parabolic cylindar functions 
with their asymptotic expansions*^ >1‘J for large 
arguments given the approximate high frequency ex- 
presotons for the far flotd pressure. The results 
are as follows. 

For the aourco in the uniform flow 


(41) 

To obtain far field pressure solutions for 
multipolo sources located In tile uhoar layer the 
equation for wave propagation in a tronavoracly 
sheared mean flow must be considered. This is 


1 Im 





l sin a 

t.</90° 

Sirrf 

0 

■ |(M 0 cos 

* • 4 


-il/2 




cod t») 2 -coB 2 i><-6+Ji)+r} +lut 


to 

1 V *i i ip, CXp 

(l - H q coa '■)((! ■ M q coo ny - cob 


16 


2H coo O 
o 


f so "l 1/2 2 r r 2 7 \ */^~ 

^(* ■ cos «)“ * COO-. J - coa w *n ul! o cob <• - 1) + ul - M q cos h) - caa*(A 


+ sin ■.* + cos In (Sin (» - 1) 


For the source in the shear layer 


(3b) 




t nln ; 

1/2 

- 1— r+i'j* 
c 

ie 0 

r *6^ 

v o 

p lm 

u/9Q u 

dirr ~ ■ 
c o_ 

(l - £ M o cos '•) [(i - | M 0 coa 

z ll/ft 
- con n 

2H con > 

L « 

r. 

h .. 

aT. 

f li \2 2 l l/2 2 

[/ h , A 

(t. h 

I 

coa " - v [( 

1 - - M q co 3('J - cos !■ - cos " In 

!( o M 0 to ° ■ 0 

+ (0 ■ S 



Note that again these results reduce to the same 
oxprcBBion when the source Is located at the edge 
of the shear layer joining the uniform flow; that 
ia, t * h. 


Hultipole Sources 


+ Bin 1 


2 

+ cou I 1 


ln(oin '* - 1)> 


(39) 


given in Chapter l of Coldstoin*4 fn a general 
form as : 



1 D 

2 2 P 

c o Dt 


This analysis can be extended to multlpolo 
aourcen located either in the uniform flow or in 
the shear layer. For o multipole line source lo- 
cated in the uniform flow region, the far field 
pressure msy be obtained by differentiating the 
monopole solution with respect to source location 
(Ref. 12). Since the source location is only de- 
pendent upon the transverse coordinate direction 
(y - direction) only lateral multipole sourcca will 


- 2 



dy ix >y 



2 


<t 

dt! * 

dy ix 



(42) 


whore I in an externally applied volume force 
arising from the momentum equation and q la an 
external volume flow source within the fluid which 
arisen from the continuity equation- If the vol- 


6 


ume flow source is put equal to taro (l.c., no mass 
la added to the flow) the two romain'fB terms on 
the right hand oido of equation (42) are due to tho 
externally applied force T, Goldstein^ shows 
that In flows with no moan voloclty an oxeornal 
'urea behaves as a volume dipole nource of strength 


Now a dipole oaurco 1 b not uniquely defined In 
a tranovoroely sheared mean flow, If we want to 
oxtond Ito definition so that Its ntrength still 
colncldeo with an external force wo con doflne It 
bo that T on tho right hand aide Of equation (42) 
la the strength of a general dipole. 

Similarly a quadrupole can be defined bo that 
It corresponds to an external b tress. This is done 
by replacing fT vLth • T where T Is now a 
tensor, Hence In particular the source tarn re- 
presenting a lateral quadrupole of atrongth IT In 
tho transverse direction 1 b 

. du <* 2 r 

Dc 2 ‘ dy HZ )y 

Hence, the solution to equation (42) with q - 0 
can be obtained by differentiating the monopole 
solution with respect to source location h once 
for the dipole (aaaoclatcd with the force vector) 
and twice for the quadrupole (oDBoclatcd with the 
stress tensor), 


Only tho more interesting quadrtpolo result 
will be given hero. 



(43) 


where p, m Is given by equation (33). This ex- 
pression gives the far field preasuro for a point 
lateral quadrupole source located In the shear 
layer. 

Inspection of equation (43) Indicates that the 
first term Inside the curly brackets will dominate 
at low frequencies while the second term (contain- 
ing ('J'/c 0 )*)) will dominate at high frequencies. 
Unto tho product of these respective terms with the 
low and high frequency approximations for p. will 
give thq low and high frequency approximations for 
the far ^leld quadrupole pressure pattern, respec- 
tively, 

’flie low frequency approximation is 
2^2 cos .) 2 

p lq (1 - (h/. )M cos ) p Im 
whore p lm is given by equation (37) 


Tho high frequency approximation is 
p tq ' M^d 0 ) 2 t(l - (h/6)M o cos O) 3 

- cos 2 0(l -(lt/6)M o cos 0>]p lm (45) 

whurc p lm Is given by equation (39). 

Tho dotulls of the proceeding analysis may bo 
found In Ref. 12, 


Results and Dlucuanlon 

Thu for field pressure pattern for sound ra- 
diating from a given type of point source in a uni- 
form flow or a ltnoor shear lnyor may bu controlled 
by a number of different factoro. Among these are 
Mach number of tho uniform flow, frequency of tha 
sound, shear layer thickness and source location 
with respect to the shear layer. The effects of 
varying these parameters are determined from tho 
appropriate pressure solutions given in tho pre- 
ceding section. It should be noted that three 
dimensional effects and turbulent scattering were 
Ignored In obtaining the pressure expressions. 

The sound pressure field was calculated by 
programming the various preasuro expressions using 
standard Fortran techniques. The parameter r has 
been normalized to 0,305 motor, The parabolic cy- 
linder functions ocurring In the exact prcooure ex- 
preoslonu were calculated using a computer program 
which was originally dovcloped In connection with 
other research boing performed at NASA Lewis Re- 
search Center. Tills program was capable of cal- 
culating up Lo 80 towns in the series for parabolic 
cylinder functions with complex arguments. At 
higher frequencies this proved to bo Inadequate. 
However, at the frequencies where computational 
difficulties were encountered with the exact oolu- 
tlons the approximate expressions had already shown 
good agreement with tho exact solutions • Theuc 
calculations wore made using single precision 
arithmetic which gave nine significant decimal 
digits. The directivity calculations were made at 
5 degree intervals. 

Figures 3 to 7 shew prcooure level referenced 
to ,>c§ (i.c., Pressure Level 10 logjo i PlA' ,c o 2 
plotted against ^ 180° - •' where Is the 
angle from the positive x-axls (downstream flow 
direction). Those curves are for a source of unit 
strength, a uniform flow Mach number of 0.8 and a 
shear layer thickness of 0.153 meter. The plots 
compare the exact and approximate calculations de- 
rived In the previous sections as indicated in the 
various figures . 

Source in the Uniform Flow 


Figure 3 shows results tor a source In the 
uniform flow. Wien treating a point source located 
In the uniform flow, the source is Considered to be 
within a few wavelengths of the shear layer. If 
the source It, many wavelengths from the shear layer 
the problem degenerates to the plane wave cose al- 
ready treated by Graham and Graham'. 

One of the most significant characteristics of 
the prcsoiiro pattern In Fig, ,i i<j the designated 
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"tana at tl lonco. 11 TUla cotton la d etermined by tlio 

quantity ^/(L - Mg cofl H) 2 - cos 2 0 which appears 
In Q*~ (Eq i OO)) as wall as In the exponent of tho 
fat field pressure oxpreaoion (Eq, (29)). The angle 
at which thin square root bocomca taro la roforred 
to as the critical ai.^ln and Is defined by 



Thlo la the srooltcst angle, o t that the transmitted 
wavo can make with the x axis. Note that this 
angle lo completely controlled by the Mach number of 
the uniform flow. Tho "tone of silence" lo defined 
as the region where <) lo leas than fl c . In this 
region the above aquaro root term acts no on expo- 
nential damping factor. Hence thu "zone of stlcnco" 
lo actually a zone of attenuated sound propagation. 

For tho region In which *> Is greater than 
'.' c the square root term contributes only to tho 
phase of tho sound wave and no longer affects the 
amplitude, 

'the far field prosauro patterns for dipole and 
quadrupole sources in tho uniform flow are shown in 
Figs. l\ and 5. Tlieae for Held pressure patterns 
also exhibit a well defined zone of silence. How- 
ever, rothor than poohlng at the critical nnglo tho 
pressure level goes to -»• due to tho fact that 
the above term to a multiplicative factor In the 
preen ire express ions (Eqa, (40) and (41), respec- 
tively) which therefore go to zero at '> c in each 
case, 

'ihe curves oliown in Figs, 3, 4, and 5 ore ail 
ff* a Urouhal number of l which means that the 
oopust.c wavelength lo approximately equal to the 
nll-'-r fiver thickness, Additional calculations 
wore made for Strouhal numbers ranging from 0.25 up 
to about 1,0. Far Strouhal numbers much above 3.0 
thu computer program for calculating the parabolic 
cylinder functions encountered computational diffi- 
culties for angles near 90°. Hence, the exact ex- 
pression for tiie sound field Is reliable only for 
certain regions of the sound field at higher iro- 
qucnciea. Outside of tho zone of silence both of 
the approximate calculations agree with the exact 
result over the range of Strouhal numbers consid- 
ered. Tills behavior indicates that neither the 
radiation pattern nor the method of calculation are 
particularly sensitive to frequensy outside of the 
zone of silence. Tills corresponds to the assertion 
of previous researchers?’ 11 that the shear layer 
thicUnoue has essentially no effect on the sound 
radiation pattern except inside the zone of silence 
( * - c ). These characteristics carry over direct- 

ly to the dipole and quadrupole calculations. 

Finally note that the pressure level tends to 
- at 0° and 180°. Thin corresponds to zero 
pressure at these anglea. This in attributed to 
the cancellation of incident woven by reflected 
waves at small incidence angles. Thin is what 
Gottlieb^ refers to as the "Lloyd 1 ;; Mirror Effect." 

Source in the Shear Layer 

The sound radiation field due to a volume flow 
source embedded in the shear layer is similar in 
nature to that of o monnpole source located in the 
uniform flow. A typical radiation pattern for a 


a 


monopolc source In tho ahoar layer la ohown In 
Fig, 6 for a Strouhal number of 1, 

For tho source In the sheaf layer there at'll 
ap pears to bo a zono of a lienee again determined by 

fit - MqCOS 0)2 - coo f 0 r tho oxact and low 
frequency resulta; however. It enters the pressure 
expression (Eq. (38)) only through tlio coefficients 
of the parabolic cyilndor functions given by G** 
(Eq. (30)) and doss not entor directly Into tho ex- 
ponent of tho oxact pressure equation or tho low 
frequency result, Hence there is no exponential 
damping factor In olthor of these equations. How- 
ever, tho combination Of algebraic terms. 

(I - M 0 cas o)2 
(l -(h/6)M 0 coa n) 2 


x nln " 

(L - H 0 coa C0 2 sin it •( f^(l -^Iptoo f)) 2 - cos 2 " 

appearing In tho low frequency pressure expression 
(Eq. (37)) results In a rapid dropoff of the sound 
levels for angles leas than ’) c « con" l 1/(1 + Mg). 
This behavior is much like the exponential , -.imping 
in tho zone of silence at high frequencies or when 
tho source is in the uniform flow, it should bo 
noted that this rather remarkable feature la In- 
depenjent of frequency while the exponential damp- 
ing factors found in the high frequency approxi- 
mation for the saute*) in the sheet layer and in all 
results for the source in the uniform flow do de- 
pend on frequency. 

The poult In the low frequency approximation 
curve occurs ot c (l.o., when tlio square root in 
the denominator of tile second factor is zero). 

As 1 decreases below -i , the denominator of the 
second factor approaches I (l.e., the square root) 
while sin . approaches zero thus causing Lhe de- 
crease In sound level as approaches zero. At 
angles near 90° sin ■ is near I causing the sound 
level to "flatten out" In thin region. 

The square root ( ^(1 - M n ioo ) 2 - cos 2 
does not appear in the high frequency approximation 
(Eq. (39)). However this approximation does con- 
tain tho square root ( f(l - (h/A)M l ,coa d*) 2 - coo 2 ") 
both in the denominator and in the exponent. Hence 
this square root determines a "critical angle" for 
the high frequency expression given by 

* c |,f “ (1 -I (lt/6)Mo) 


(hose that (it/e)!J 0 Is the local Mach number of the 
flow at the source locstfon in the shear layer.) 
this accounts f. r the fact that tin peaks fur the 
high frequency approximation occur at slightly dif- 
ferent angles than the peaks ul the exact and low 
frequency curves- 

IL should he noted that both of the square root 
e xpressions arise o riginally tram the expression 

Vk 2 - ( /c 0 - Ml) 2 i( /c„)V(l - H coa - )2 - coo 2 **. 
Hence the square* rout Is zero at branch points of 
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the complex k-pisno and as shown above these points 
determine "critical angles." It Is at these points 
that the high frequency expansion brooks down. 

Ti e idoo of a dipolo source or a quedrupolo 
source (or for that matter any ocher order multi- 
polo source) In a transversely sheared mean flow Is 
rather vague since thoro is no unique way of de- 
fining such a source. In the previous section 
thcac sources ore defined so that they correspond 
to an externally applied force and on externally 
applied stress, respectively. The curves shown in 
FIs. 7 are for tha quodrupola source located at the 
center of the sliaar layor, a neon flow Mach number 
of 0,S and Strouhal numbers of 0.5, 1.0, and 6.0, 
rospoctivoly. 

The exponential damping is most easily seen In 
tho high frequency approximation glvon by aquation 
(39). However, this feature Is not as pronounced 
as might be expected. Tills is because of a masking 
resulting from tho coefficient In equation (65). 

One further point of significance which applies 
In general to this typo of quodrupolo is that the 
high frequency approximation equation (66) resembles 
tho exact express ion for a quodrupolc In tho uni- 
farm flow equation (65) multiplied by n dopplur 
factor (1 - h/ti HqCos <>), where tho Mach number In 
each esse is that of the flow at tho source loca- 
tion. 


Summary and Conclusion 

An analysis has been developed for calculating 
the far field acoustic pressure for sound waves pro- 
pagating through o linear shear layer. Tho sound 
waves arc produced by a line source either located 
nosr the shear layer in a uniform flow or located 
in tho shear layor Itself. Closed form analytic 
express Iona for the far field pressure have been 
obtained for several different types of sources 
(l.o. , monopole, dipolo, etc.), Tho definitions of 
these sources depend upon whether the source is lo- 
cated in the uniform flow or in the shear layer. 
Approximate ei ivcssiens are given for limiting cases 
of high and low frequencies for each type of source. 
These approximate expreasions give very good agree- 
ment with the exact result tor cheir respective 
limiting cases. 

The exact pressure expression for a monopole 
source reduces to previously known results for such 
special esses ns zero shear layer thickness, pro- 
pagation of the sound wavog perpendicular to the 
flow direction, and plane waves propagating through 
a velocity discontinuity, The exact expressions 
for multipole sources in the uniform flow are sim- 
ple extensions of the monopole result. The expres- 
sion for tha monopole or volume flow source In the 
shear layer becomes Identical with the monopole ex- 
pression in the uniform flow when tho source la 
placed at the edge of the shear layer. The quadru- 
ple Bourcc solution in tho shear layer is an ex- 
tension of the monopole solution in tho shear layer. 
Consequently, all of the pressure expressions de- 
veloped hero can be reduced to known results lor 
special limiting cases thus providing a check for 
the analysis. Hence analytic expressions have been 
developed far calculating the sound field due to 
various typcB of line sources located in a uniform 
flow or in a shear layer for the entire acoustic 


frequency spectrum. The asymptotic expansions of 
these exact expressions for tho limiting cases of 
low and high frequency correspond to results pre- 
sented elsewhere for sheer leyers of arbitrary 
velocity profiles. 
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Figure 1. - Propagation of aircraft nolso through shear layers. 
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Figure 2 . - Schematic of simplified shear layer used in analysis. 




PRESSURE LEVEL -10 


c\i 

Cvj O 
u 



«— I 

S’ 


EXACT CALCULATION {EQ. {29)) 

HIGH FREQUENCY APPROXIMATION (EQ. (33)) 

LOW FREQUENCY APPROXIMATION (EQ. (36)) 



SILENCE 


-145- 1 1 ■■ 1 - 1 ; 1 i 1 i 

0 20 40 60 80 100 120 140 160 180 


0' = 180° - 0 


Figure 3. - Comparison of exact and approximate 
pressure expressions for a monopole source 
in the uniform flow. 


EXACT CALCULATION (EQS. (40)) 

HIGH FREQUENCY APPROXIMATION (EQS. (40); 

LOW FREQUENCY APPROXIMATION (EQS. (40)) 



Rgure 4. - Comparison of exact and approximate 
pressure expressions for a transverse dipole 
source in the uniform flow. 


■ EXACT CALCULATION (EQ. {41]) 

HIGH FREQUENCY APPROXIMATION (EQS. (41)) 
LOW FREQUENCY APPROXIMATION (EQS. (41)) 



8' = 180° - 0 



Figure 5. - Comparison of exact and approximate 
pressure expressions for a transverse quad- 
rupole source in the uniform flow. 
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EXACT CALCUiATlON (EQ. (33) 

HIGH FREQUENCY APPROXIMATION (EQ. (39)) 



8' = 180P - 8 

Figure 6. - Comparison of exact and approximate 
pressure expressions for a monopole source 
in the shear layer. 




(ACT CALCULATION (EQ. (43)} 

GH FREQUENCY APPROXIMATION (EQ. (44)) 
>W FREQUENCY APPROXIMATION (EQ. (45)) 





